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Symmetry potential of A(1232) resonance and its effects on the vr /7r+ ratio in 
heavy-ion collisions near the pion production threshold 
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Effects of the completely unknown symmetry (isovector) potential of the A(1232) resonance on 
the total and differential 7r“/7r’^ ratio in heavy-ion collisions at beam energies from 100 to 1000 
MeV/A are explored within an isospin-dependent transport model IBUU. The effects are found to 
be negligible at beam energies above the pion production threshold due to the very short lifetimes of 
less than 2 fm/c for A resonances with masses around mA = 1232 MeV, leaving the ratios of 

especially the energetic pions still a reliable probe of the high-density behavior of nuclear symmetry 
energy Esym{p)- However, as the beam energy becomes deeply sub-threshold for pion production, 
effects of the A symmetry potential becomes appreciable especially on the n~ ratio of low-energy 
pions from the decays of low-mass A resonances which have lived long enough to be affected by their 
mean-field potentials, providing a useful tool to study the symmetry potential and spectroscopy of 
A resonances in neutron-rich nuclear matter. Interestingly though, even at the deeply sub-threshold 
beam energies, the differential 'k~ ratio of energetic pions remains sensitive to the Esym{p) at 
supra saturation densities with little influence from the uncertain symmetry potential of the A 
resonance. 

PACS numbers: 21.65.Ef, 24.10.Ht, 21.65.Cd 


I. INTRODUCTION 

Thanks to the great efforts of many people in both 
nuclear physics and astrophysics communities, much 
progress has been made in recent years in constrain¬ 
ing the density dependence of nuclear symmetry en¬ 
ergy E^ymip) around the saturation density pa, see, e.g., 
refs. [iHllI for comprehensive reviews. At supra sat¬ 
uration densities, however, even the tendency of the 
Esym{p) is still controversial based on the model analy¬ 
ses of very limited data available. Moreover, predictions 
on the high-density behavior of nuclear symmetry energy 
based on both microscopic many-body theories and phe¬ 
nomenological models diverge broadly mainly because 
of our poor knowledge about the spin-isospin depen¬ 
dence of nuclear many-body forces, isospin-dependence 
of short-range correlations, the nature of tensor forces, 
etc, besides the well-known difficulties of treating pre¬ 
cisely quantum many-body problems m- 

Heavy-ion collisions are the only means in terrestrial 
laboratories to create dense nuclear matter. To pin down 
the high-density symmetry energy has been one of the 
major goals of the low-intermediate energy heavy-ion re¬ 
action community in recent years. Among the potential 
probes of the high-density symmetry energy, the 'k~ j'K'^ 
ratio in heavy-ion collisions has attracted much atten¬ 
tion since it was first propo sed (T^ . Unfortunately, sev¬ 
eral model analyses [l3l - ll8l | of some existing data have 
so far been inconclusive. Realizing the importance of 
thoroughly understanding all aspects of the tt — N — A 
dynamics in isospin-asymmetric matter especially since 
several experiments measuring the 'k~ ratio are cur¬ 
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rently underway at RIKEN [1^ and NIRS [1^ besides 
those planned at other facilities [lil, significant ef¬ 
forts have been made recently in studyiM effects of the 
in-medium pion production threshold |23j . pion mean- 
field [ 234 ?^ . electromagnetic field [ 2 ^, neutron-skin [ 2 ^ 
and nucleon-nucleon short-range correlation HEii on 
pion production mostly at intermediate energies above 
the pion production threshold of about 300 MeV/A. Ex¬ 
tending these efforts, we report here results of a study 
on effects of the completely unknown symmetry (isovec¬ 
tor) potential of the A(1232) resonance on both the to¬ 
tal and differential tt ~ratios in heavy-ion collisions 
at beam energies from 100 to 1000 MeV/A within an 
isospin-dependent transport model IBUU [l^ • It is found 
that the 7r“/7r+ ratio of energetic pions remain a useful 
probe of the high-density Esym{p) without much influ¬ 
ence by the completely unknown symmetry potential of 
the A resonance. 


II. MODELING THE SYMMETRY 
POTENTIALS OF NUCLEONS AND A(1232) 
RESONANCES 

We first recall here briefly the major ingredients and 
new modifications of the IBUU transport model that are 
relevant for this study. First of all, to ease future compar¬ 
isons with studies using other transport models, for the 
isoscalar potential we use the momentum-independent 
BKD (Bertsch-Kruse-Das Gupta) nucleon potentials 
with an incompressibility of ko = 236 MeV that is a 
common option in most transport models developed for 
investigating heavy-ion reactions. This choice does not 
affect our qualitative conclusions as we are dealing with 
the ratios of charged pions that are insensitive to the nu¬ 
cleon isoscalar potential. On the other hand, it allows us 
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FIG. 1: (Color online) The density dependence of the poten¬ 
tial symmetry energy used in this study. 


to generate timely enough reaction events necessary for 
studying various cases of deeply sub-threshold pion pro¬ 
duction using the best high-performance computing re¬ 
sources available to the author today. We caution, how¬ 
ever, to compare with data of ongoing experiments to 
constrain quantitatively the Esymip) at supra saturation 
densities using the differential ■n~ j-K^ ratio of energetic 
pions, the momentum dependence of both the isoscalar 
and isovector potentials of nucleons and A resonances 
should be considered carefully. 

For the purpose of this work, it is sufficient and more 
economic to use the simplest parameterization for the 
density dependence of nuclear symmetry energy Esym (p) ■ 
Since the earlier studies about effects of nuclear symme¬ 
try energy on several properties of neutron stars ISp . 
and the isospin dynamics of heavy-ion collisions |35l | , one 
often parameterizes the Esym{p) as 

E,ym{p) = (2^/3 _ [«2/3 _ ^ E,y^{po)F{u)il) 

with the typical function F{u) being 

27/2 

= T—> -^2(w) = u, Fsiu) = (2) 

1 -I- u 

where u = p/po is the reduced baryon density, Ep = 36 
MeV is the Fermi energy and Esymipo) = 31 MeV is the 
symmetry energy at po. The selected three forms of the 
F(u) span a large uncertain range of Egymip)- What is 
important for the present study is the potential symme¬ 
try energy given in the second part of Eq. [T] and shown 
in Fig. [TJ Assuming there is no momentum dependence, 
the nucleon symmetry potential Vg^syin/p) up to the lin¬ 
ear term in the isospin asymmetry S = [pn — Pp)/iPn+Pp) 


of nucleonic matter can be obtained from taking the par¬ 
tial derivative of the potential symmetry energy density 
with respect to the neutron or proton density The 
Va.sy(n/p) corresponding to the three functions of F{u) is 
given by, respectively, 

Vasyl(n/p) = idea—- S, 

1 + u 

Vasy 2 in/p) = ±2eaUS, 

Vasy3{n/p) = ±2eaU^^‘^S (3) 

where the ± sign is for neutrons/protons and Ca = 
Esymipo) — (2^^^ — 1) ^Ep. We notice here that the po¬ 
tential symmetry energy density for Fi(u) and F^^u) also 
contributes a term proportional to S'^ to the isoscarlar po¬ 
tential. However, such a term in heavy-ion collisions is 
negligibly small compared to the isoscalar potential cor¬ 
responding to the EOS of symmetric matter. 

Associated with the tt-N-A dynamics in nuclei, nu¬ 
clear reactions and neutron stars, the isoscalar potential 
for the A(1232) resonance has been extensively studied 
for a long time using various many-body theories and 
interactions, see, e.g., refs. [36l - l4^ . For a comprehen¬ 
sive review of the relevant issues, we refer the reader to 
ref. [i^. Indeed, there are indications from analyses of 
electron-nucleus, photoabsorption and pion-nucleus scat¬ 
tering that the A isoscalar potential va is in the range 
of —30 MeV + un < UA < un with respect to the nu¬ 
cleon isoscalar potential un [d^ . Since a constant poten¬ 
tial has no dynamical effect in heavy-ion reactions and 
there is no information available as to how the differ¬ 
ence in the isoscalar potentials for nucleons and A reso¬ 
nances depend on the density and/or other properties of 
the medium, in this study we thus use the same isoscalar 
potentials for nucleons and A resonances. The mean- 
field potential of pions in neutron-rich matter is an in¬ 
teresting issue. While there is some information about 
the strength of pion potential in nuclei from study pio- 
nic atoms, there is little information about its density, 
isospin and momentum dependence in nuclear matter. 
Recognizing its possible importance but also considering 
the currently existing controversies about it, we neglect 
here the pion mean-held. 

As pointed out already in ref. [3 IH, 13 > essen¬ 
tially know nothing about the isovector (symmetry) po¬ 
tential of the A(1232) resonance in isospin-asymmetric 
nuclear matter. In heavy-ion collisions near the pion 
production threshold, the A population is very small. 
For example, for ^^^Sn-|-^^^Sn reaction at a beam energy 
of dOO MeV/nucleon and an impact parameter of 1 fm, 
the maximum total multiplicity of A resonances in all 
charge states reached during the reaction is about d 3 ■ 
This is about 3.2% of the total number of nucleons in¬ 
volved. Thus, the mean-held potential of A is mainly 
due to its interactions with nucleons not other A reso¬ 
nances. Therefore, similar to nucleons, the isovector po¬ 
tential of A resonances is also proportional to the isospin 
asymmetry 6 of nucleons due to the isospin-dependent 
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T 3 (A)-r 3 (iV) term in the A—N interaction [i^. The pop¬ 
ulation of A only has a small perturbative effect on the 
isospin asymmetry S of nucleonic matter. Assuming the p 
meson exchange among nucleons and baryon resonances 
is mainly responsible for the isovector interactions, there 
is a simple scaling relation between the NNp and AAp 
coupling constants (?pN and q„/\ in tree-space based on the 
quark model of baryons [3g. However, the relationship 
between the isovector potentials of nucleons and baryon 
resonances in isospin-asymmetric matter is completely 
unknown. In fact, the critical density for A(1232) res¬ 
onance formation in neutron stars has been found to de¬ 
pend almost linearly on the unknown ratio Xp = PpA/PpN 
[^1^. Depending on the assumed value of Xp, the for¬ 
mation of A(1232) resonance can even happen a t pn and 
affect significantly properties of neutron stars |46l |. In 
addition, the isovecor Lorentz-scalar i5-meson may also 
play some roles in determining the baryon isovector po¬ 
tentials. However, its couplings with baryons are even 
less known. Thus, because of their ramifications in both 
nuclear physics and astrophysics both the isoscalar and 
isovector parts of the A potential deserve further inves¬ 
tigations. In particular, it would be interesting to know 
if any observable in heavy-ion collisions especially those 
induced by neutron-rich nuclei may be used to probe the 
A isovector potential in neutron-rich matter. 

In two earlier studies [mill, the isovector potential 
of A resonances was connected to that of nucleons based 
on rather different assumptions. In ref. [H and sub¬ 
sequently many other studies, considering the A as a 
molecule consisting of a nucleon and a pion, the isovec¬ 
tor potential of the A resonance is an average of that for 
neutrons and protons with weights given by the square of 
the Clebsch-Gordon coefficients in the A o irN processes 
conservii^the total isospin. Under these assumptions, 
one has [1^ 

'Casy(A ) = Vasy{j^\ 

2 11 

Vasy{A°) = -Vasy{n) + -Vasy{p) = -^Vasy{n), 

12 1 

'^asy{A ) — T ~^'^asy{p') — 


Vasy{A++) 


Vasy{p) = -Vasy{n). 


(4) 


Assuming that baryon isovector potentials are all due to 
the p meson exchange neglecting possible contributions of 
the (5-meson, the above prescription is equivalent to set¬ 
ting Xp = 1/3. In line with the above assumption, taking 
into account the small contributions of the A resonances 
the effective isospin asymmetry Suke of the excited bary- 
onic matter can be obtained from 


^like — 


{Pn'jlike (Pp) like 
(,Pn')like + (pp) like 


where 

/ N 2 1 

[Pnjlike — Pn + gPA^ + 3PA+ + PA” , 

( N 2 1 

[PpJUke — Pp T 2PA+ 4 “ 3^A0 T Pa++- 


(5) 

( 6 ) 
(7) 


It is evident that the Suke reduces naturally to the isospin 
asymmetry S of nucleonic matter as the beam energy 
becomes smaller than the pion production threshold. 

While in ref. , the following prescription according 
to the effective charge of the A(1232) resonance was used 


'^asy{A ) 
Vasy{A°) 
‘Casy (A^) 

Ua.y(A++) 


2Vasyin) - Vasyip) = 3Vasy{n), 
Vasyin), 

Vasyip) = -Vasyin), 

‘^Vasyip) Vasyi^l) — ^Vasyi^l) • (8) 


This prescription is equivalent to setting Xp = 1, thus 
the magnitudes of these isovector potentials for A(1232) 
resonance are 3 times that of Eq. SI We emphasize that 
the ratios of A isovector potentials for different charge 
states are the same in both prescriptions. In this work, 
we use the symmetry potentials in Eq. SI as a basis and 
multiply them with a “A-probing factor /a” . By varying 
the value of /a we can explore effects of the isovector 
potential of A(1232) resonance on pion observables in 
heavy-ion collisions. 
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FIG. 2: (Color online) The mean lifetime (black line) of 
A(1232) resonance as a function of its mass. The dashed line 
at 1.7 fm/c indicates the lifetime of the A as its centroid mass 
of 1232 MeV. The blue lines are the A mass distributions in 
NN —>■ NA collisions at nucleon-nucleon (NN) center of mass 
energy of =2.11 GeV and 2.32 GeV corresponding to first 
chance NN collisions in nucleus-nucleus reactions at a beam 
energy of 500 and 1000 MeV/nucleon, respectively. 







III. MODELING THE MASS DISTRIBUTION 
AND LIFETIME OF A(1232) RESONANCES 
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Whether the A potential has any observable effect in 
heavy-ion collisions depends strongly on the lifetime of 
A(I232) resonance over which the A potential acts to 
modify its momentum. It is thus important to outline 
how we model the mass distribution and the correspond¬ 
ing width of the A(I232) resonance in heavy-ion colli¬ 
sions. As discussed in detail in ref. and reviewed in 
ref. [ 5 ^, in the IBUU transport model the mass of a 
baryon resonance r (A or N*) produced in an inelastic 
nucleon-nucleon collision {N + N ^ N' + r) is gene rated 
according to a modified Breit-Wigner function 541 


P{mr) 


Pf ■ rrir ■ T{mr) 


(9) 


where p/ is the final nucleon momentum in the NN center 
of mass frame, m^o and r(mj.) are the centroid and width 
of the resonance, respectively. As noticed earlier , the 
pf surpreses significantly the A production near its max¬ 
imum mass of mA(niax) = ^/s — ttin compared to the 
simple Breit-Wigner function widely used in many trans¬ 
port models. This effect is especially significant when 
the mA{max) is less than 1232 MeV in sub-threshold 
pion production. For r = A(I232), the width r{mr) is 
given by 

r(A(I232)) = 0.47gV {^1 + O.bg^) (GeV) (10) 


where q is the pion momentum in the A rest frame in the 
A ^ TT + N decay process. During each time step of very 
small length dt, a Monte Carlo sampling of the A decay 
is carried out according to the probability 


Pdecay = I - exp{-dt/TA) dt/TA (H) 

where ta = h/T{mA) is the A mean lifetime. 

It is necessary to mention that there are many interest¬ 
ing issues regarding the A spectroscopy in neutron-rich 
matter to be studied [i^, such as the possible modifi¬ 
cation of its width (lifetime) and centroid mass. More¬ 
over, there are already some interesting work on the de¬ 
lay times oi NN ^ NA scattering, how the width of 
A(1232) resonance in the scattering should be modeled 
quantum mechanically and whether pions should be di¬ 
rectly produced without going through the intermediate 
A(1232) resonance at all in heavy-ion collisions [sgI - Is^ . 
These issues are not addressed here. Our calculation of 
the lifetime of the A(1232) resonance is rather standard 
used in most transport models. Moreover, as discussed 
in detail in refs. [4^ [50j| . in IBUU we use free-space 
NN ^ NA cross sections and the scattering kinematics 
is that in vacuum. This is a good approximation when 
nucleons and A resonances have the same momentum- 
independent potential, since in this case the baryon po¬ 
tentials before and after the NN collision cancel out in 
the energy conservation equation. In other cases, how¬ 
ever, our treatment is equivalent to neglecting residual ef¬ 
fects of the baryon mean-field potential on the NN center 
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FIG. 3: (Color online) The {tt~/ n^)iike ratio as a function of 
time in head-on An-|-Au reactions at a beam energy of 100, 
200, 500 and 1000 MeV/nucleon with the same F(u)=u but 
3 different /a values. 


of mass energy ^/s, except those through the momenta 
of the colliding nucleons from their earlier propagation 
in the medium. As we noticed earlier, how the isospin 
and momentum dependent potentials of nucleons and A 
resonances may affect particle production thresholds is 
an interesting issue under investigation by several other 
groups. We also plan to address this issue within the 
IBUU approach in the near future. 

As discussed earlier and seen clearly in Fig. [21 the 
low-mass A resonances predominantly produced in NN 
collisions near the pion production threshold have very 
long lifetimes during which the A mean-field may act. 
Pions from decays of these resonances may then carry 
some information about the A potential. On the con- 
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kinetic energy in head-on Au-|-Au reactions at a beam energy 
of 200 MeV/nucleon with the same F(u)=u but three different 
isovector potentials for the A resonance. 
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FIG. 5: (Color online) The (7r“/7r^)iifee ratio as a function of 
time in head-on Au-I-Au reactions at a beam energy of 200 
MeV/nucleon with /a — I but two different forms for the 
F(u). 


trary, as indicated by the dashed line in Fig. [2l the A 
resonance around its mass centroid of 1232 MeV has a 
lifetime of only about 1.7 fm/c which is too short for 
the A mean-held to cause any signihcant change in its 
momentum. Moreover, subsequent tt + JV ^ A reaction 
cycles tend to destroy the memory of the pion about the 
dynamical history of its parent A. The two blue lines are 
the A mass distributions in NN —)■ NA collisions at the 
NN center of mass energy of i/s = 2.11 GeV and 2.32 
GeV corresponding to the hrst-chance NN collisions in 
nucleus-nucleus reactions at a beam energy of 500 and 
1000 MeV/nucleon without considering contributions of 
the initial nucleon Fermi momentum, respectively. It is 
seen that at 1000 MeV/nucleon, the A mass distribution 
peaks around 1232 MeV. The A resonances created at 
this beam energy are thus mostly very short-lived. As 
the beam energy decreases towards sub-threshold ener¬ 
gies, e.g., at 500 MV/nucleon, more long-lived low-mass 
A resonances are produced. One thus expects that pions 
from heavy-ion collisions at sub-threshold beam energies 
to be more sensitive to the A potential than reactions 
at higher beam energies. Moreover, only a small fraction 
of nucleons go through the A resonances for a relatively 
short time interval compared to the whole reaction pe¬ 
riod over which the nucleon potential acts continuously. 
We thus do not expect to see any effect of the A potential 
on nucleon observables. 



(MeV) 

FIG. 6: (Color online) The {'it~/ n^)iike ratio as a function of 
pion kinetic energy in head-on Au-|-Au reactions at a beam 
energy of 200 MeV/nucleon with /a = 1 but two different 
forms for the F(u). 
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IV. EFFECTS OF THE SYMMETRY 
POTENTIALS OF A(1232) AND NUCLEONS ON 
THE TOTAL AND DIFFERENTIAL 7r"/7r+ 
RATIOS 


We first examine in Fig. [3] the evolution of the 
(7r“/7r+);ifee ratio 


(tT 


7r-+A- + |A° 
7 r+ + A++ + iA+ 


( 12 ) 


for head-on Au-|-Au reactions at beam energies of 100, 
200, 500 and 1000 MeV/nucleon with F{u) = u and 
/a = 1, 3 and 6, respectively. Approximately 120,0000 
events were generated in each case. The {'!t~ lTT'^)iike ra¬ 
tio naturally becomes the final tt~I ratio at the freeze- 
out when the reaction time t is much longer than the 
lifetime of the A resonance ta . The extremely long-lived 
and very-light A resonances with ttia near rriN +'mT^ are 
forced to decay at the freeze-out time if they were not 
reabsorbed by the NA —>■ NN reactions earlier. With 
/a = 1 and 3, the symmetry potential for the A(1232) 
resonance is that of Eq. 0] and Eq. [51 respectively. We 
also use a large value of /a = 6 to explore its effects. 
It is seen that effects of the A symmetry potential only 
become appreciable when the beam energy approaches 
the pion production threshold. At deeply sub-threshold 
beam energies, the effect of the A symmetry potential on 
the total 7r“/7r+ ratio is significant. At A^eam/A = 200 
MeV, for example, the total 7r“/7r+ ratio increases by 
about 25% when the /a is increased from 1 to 6. As one 
expects, the final {tt~/T r~^)iike ratio is higher with a larger 
/a value which amplifies the difference in potentials for 
A~ and A++. The beam energy dependence of the ef¬ 
fect of the A symmetry potential is consistent with our 
expectations based on the A lifetime shown in Fig. [5] and 
discussed earlier. Namely, only pions from the decays of 
light but long-lived A resonances created in low-energy 
collisions are affected by the A potential. Also, consis¬ 
tent with the known systematics of pion production , 
the TT~ /Tr’*' ratio increases with decreasing beam energy. 

The results shown in Fig. jSjseem to make the current 
efforts of using the sub-threshold 7r“/7r+ ratio to con¬ 
strain the Egymip) around 1 — 2po even more challeng¬ 
ing unless one knows the exact relationship between the 
isovector potentials of nucleons and A resonances. To as¬ 
sess the situation more clearly, we examine in Fig. 0|the 
final TT~ ratio as a function of pion kinetic energy for 
the head-on Au-|-Au reactions at E^eam/A = 200 MeV. 
It is interesting to see that effects of the A symmetry 
potential from varying the /a parameter is mainly on 
low energy pions below about 70 Mev. The more ener¬ 
getic pions are not much affected at all by the A potential 
even at such a deeply sub-threshold beam energy. Again, 
this is due to the mass-dependent lifetime of the A res¬ 
onance. The low (high) energy pions are mainly from 
the long- (short-) lived, low (high) mass A resonances. 
Thus, to investigate the isovector potential of A(1232) 
resonance one should focus on the 7r“/7r+ ratio of low 


energy pions for which the Coulomb fields are also im¬ 
portant but are well understood. 

The observation that the uncertain isovector poten¬ 
tial of the A(I232) resonance has almost no effect on the 
differential 7r“ /tt’’' ratio of energetic pions is very encour¬ 
aging. It is thus still possible to use these energetic pions 
to probe the Egymip) at supra saturation densities. To 
be more quantitative and make a comparison, shown in 
Fig. [5] are the /'!r'^)iike ratio as a function of time in 
the head-on Au-|-Au reactions at a beam energy of 200 
MeV/nucleon with /a = 1 but two different forms for 
the F(u). As known before, the total tt~ /tt~^ ratio at the 
freeze-out has an appreciable sensitive to the variation 
of F(u). We notice that the current uncertainty of the 
potential symmetry energy at supra saturation densities 
going from large negative to positive values is much larger 
than the difference between the Fi(tt) = 2tt^/(l -l-tt) and 
E 3 {u) = used as examples here More interest¬ 
ingly, as shown in Fig. [5] the differential ratio of 

energetic pions still shows an appreciable sensitivity to 
the variation of F(u). Since these energetic pions are 
not much influenced by the unkonwn A potential, the 
observations here similar to the findings from studying 
the differential 'k~ j'K'^ ratio in Sn-|-Sn reactions at simi¬ 
lar beam energies within the PBUU model [l^, indicate 
that the 7r“/7r+ ratio of energetic pions is still a useful 
probe of the Egymip) at supra saturation densities. 


V. SUMMARY 

In summary, within the IBUU transport model we 
studied effects of the completely unknown A symme¬ 
try potential on the total and differential 7r“/7r+ ra¬ 
tio in heavy-ion collisions at beam energies from 100 to 
1000 MeV/A. We found that the A symmetry poten¬ 
tial has negligible (significant) effect on the integral ra¬ 
tio of charged pion multiplicities at beam energies above 
(below) the pion production threshold. In sub-threshold 
reactions, the differential ratio at low kinetic en¬ 

ergies is also significantly affected by the A symmetry 
potential, providing a useful tool to study the symmetry 
potential and spectroscopy of A resonances in neutron- 
rich nuclear matter. On the other hand, the differential 
7 r“/7r+ ratio of energetic pions is little affected by the 
A symmetry potential even at the deeply sub-threshold 
beam energies, leaving it still one of the best known 
probes of the Egymip) at supra saturation densities. We 
also found that the mass-dependent A width (lifetime), 
i.e., the A spectroscopy in neutron-rich matter, is very 
important in determining whether the A potential plays 
any role in heavy-ion collisions. 
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